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Abstract: In this paper we calculate the transport coefficients of 
plasmas formed by air and water vapor mixtures. The 
calculation, which assume local thermodynamic equilibrium 
(LTE) are performed in the temperature range from 500 to 
12000 K. We use the Gibbs free energy minimization method to 
determine the equilibrium composition of the plasmas, which is 
necessary to calculate the transport coefficients. We use the 
Chapman-Enskog method to calculate the transport 
coefficients. The results are presented and discussed according 
to the rate of water vapor. The results of the total thermal 
conductivity and electrical conductivity show in particular that 
the increasing of the rate of water vapor in air can be 
interesting for power cut. This could be improve the 
performance of plasma during current breaking in air 
contaminate by the water vapor. 

Keywords: Plasma, Equilibrium composition, Transport 

coefficients. 

1) Introduction 

The current cut in the circuit breakers is obtained by 
separating electrodes in air or in compressed air [1, 2]. When a 
short-circuit or an overload appears in the electrical circuit, the 
contacts of the circuit breaker separate and generates an electric 
arc [1-3] which interacts with the surrounding gas and creates a 
plasma of this gas [1, 2]. Other studies on the electrical arc had 
shown the influence of the insulator vapor [1, 5-8], metal of 
contacts [1, 9-11] and hydrogen atom [1, 2, 5, 8] on the 
engendered plasmas characteristics. In this paper, we are 
interested by the influence of hydrogen on plasma characteristics 
for electrical current breaking. This hydrogen can be brought by 
water vapor in the air. The knowledge of transport coefficients is 
fundamental to evaluate and compare the performance of plasmas 
for current breaking [5] and also, transport coefficients are the 
indispensable data for modeling of electric plasmas arcs [12- 
14]. We undertake in this study the calculation of the transport 
coefficients of air-water vapor mixtures thermal plasmas. This 
mixture can be naturally formed because of air humidity [1, 2], 
but the percentage of water vapor in the atmospheric air is 
limited. In the zones with great humidity, we can have a 
maximum percentage of water vapor in air around 7% at a 
temperature of 40°C [1, 2]. As this value is not sufficient to bring 
enough hydrogen for better performances, we choose to study air- 
water vapor mixtures with often high water vapor percentage. The 
plasmas considered are: 95% air - 5% water vapor; 90% air - 10% 
water vapor, 50% air - 50% water vapor and 20% air - 80% water 
vapor. These four mixtures are chosen to determine the water 
vapor influence on the plasmas transport coefficients. The 



increase of water vapor in the mixture could be done during the 
breaking of the electric current by water injection. Assuming local 
thermodynamic equilibrium (LTE), calculations are made in a 
temperatures range from 500 K to 12000 K. We suppose that the 
air is constituted of 80% nitrogen (N 2 ) and 20% oxygen (0 2 ). In 
the figures 1, 2 and 3, we compare our results of dry air plasma 
with those of P. Fauchais and al. Then, we present and discuss our 
results of transport coefficients of plasma formed of air-water 
vapor mixture. This study completes the one of Kagone and al [1] 
to have a database on the range of temperatures from 500 K to 
30000 K for the modeling of plasma formed by air- water vapor 
mixtures. 

2) Transport coefficients 

The initial step in the calculation of transport coefficients is 
the determination of the equilibrium composition of the plasma. 
This is done by using the widely known method based on the 
Gibbs free energy minimization [12, 14, 16]. For air-water vapor 
mixture, we take into account 11 monatomic species: H, H, H + , 
N, N + ,N ++ ,0, O", 0 + , 0 ++ , and electrons; 18 diatomic species: H 2 , 
Hz, H 2 + , N 2 , N 2 , N 2 + , NH, NH , NH + , NO, NO , NO + , 0 2 , 0 2 , O z 
+ , OH, OH", OH + , and 1 1 polyatomic species: 0 3 , H z O, N z O, NH 3 , 
N 2 0 3 , HN0 2 , H0 2 , N 2 0 4 , H 2 N 2 , N 2 0 5 , HN0 3 . The monatomic 
and diatomic specific thermodynamic properties such as 
enthalpies, chemical potentials are calculated by the fitting 
coefficients found in [15]. For polyatomic chemical species their 
specific thermodynamic properties are obtained in Janaf tables 
[17]. 

The transport coefficients were calculated using the Chapman- 
Enskog method based on sonine polynomials development [1,5, 
12, 13, 18-22]. In our study we considered thermal conductivity, 
electrical conductivity and dynamical viscosity of the plasma. 
2.1) Thermal conductivity coefficient 

The total thermal conductivity X tot can be separated into four terms 
with a good accuracy [1,14]: 

Aot = Kr + Kr + Ant + Keact W 

Where A.* r is the translational thermal conductivity of electrons, 

A* the translational thermal conductivity of heavy species 

particles, A int the internal thermal conductivity and l react the 

chemical reaction thermal conductivity. Thermal conductivity of 
electrons at the third approximation order is given by [1, 14]: 

(2) 



InRT 
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The translational thermal conductivity due to the heavy species in 
the second order approximation can be written as [23]: 
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The coefficients Ly are given by Muckenfuss [13, 23]. 
The internal conductivity due to the effect of internal degrees of 
freedom is taken into account with the Eucken correction [14]. 
The first order in approximation of Chapman-Enskog's method 



[13] given by: 
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Where R is the perfect gas constant, Cpi the specific heat capacity 
of i th chemical specie, T the temperature and Q' im the internal 

partition function of i th specie. The formulation of the thermal 
conductivity chemical reaction was developed by Butler and 
Brokaw [13, 14, 18] and is written as: 



^react 
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i"' reaction can be written as: 



X k X, 



(8) 
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Where B is the symbol of the k chemical species. The |i 
reactions that we have to take into account must be linearly 
independent. For the i th reaction the enthalpy variation L \H i is 

given by: ^ H ^r- . „ (10) 



Where Hj is the specific enthalpy of the j chemical species and P 
the pressure. 

2.2) Electrical conductivity coefficient 

The electrical conductivity a at thermal equilibrium is obtained 
using the third-order approximation of the Chapman-Enskog 



method [14, 18]: 
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Where e is the electronic charge, ni the electrons density number, 
k the Boltzmann constant, mi mass of an electron and T is the 



temperature. The parameters q 1J are a function of density and 
collisions integral are given in [14]. 
2.3) Coefficient of viscosity 

The viscosity can be considered mainly due to the heavy chemical 
species because of the mass ratio between electrons and heavy 
chemical species. We use the formulation given in [13]: 
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Where Hy are function of collision integrals, molar mass and 
molar fraction. 

3) Choice of the potential and collision integrals 
3.1) Neutral - neutral interactions 

Many type of potential have been suggested to describe the 
interaction between two neutral particles [12]. In this paper we 
choose to use (12 - 6) Lennard-Jones potential for interaction 
between two atoms and deduce from simple combining rules the 
Lennard-Jones potentials for the other interaction for other neutral 
species. We have used the values found in the works of Andre 
and al. [14]. In table 1, we summarize the values of the Lennard- 
Jones parameters used in the case of the interaction between two 
neutral atoms simple combining rules [12] .When the parameters 
of (12-6) Lennard-Jones potentials are determined, the average 
cross sections are determined from the table of Klein and Smith 
[12, 14, 18,24]. 

Table 1: Parameters of Lennard -Jones potentials used in the 
present work. 



Species 


(°) (A) 


(£) (K/particle) 


References 


H H 


2.4 


40 


[12] 


O O 


3.050 


106.700 


[111119] 


N N 


3.298 


71.400 


[11] [12] 



3.2) Electrons - neutral interactions 

For electrons-neutral collision, the average cross sections q { -- s) 

with s included between 1 and 5, are calculated by [14]: 

(13) 



Q)j = a, ln(r) + a 2 + a,T + y T + a 5 T~ 2 + a^T' 1 + a 7 T~ 

T is the temperature, the fitting coefficients a t are taken in [14]. 
For missing data we use the rigid spheres approximations [24]. 
3.3) Ion - neutral interactions 

For M —M + interactions, the resonant charge exchange is a 
dominant process to determine the average collision integrals 
[11,14,18,23,25]: 



2r j ;r ^ 1 
4 6 



T 

2 I M 



(14) 



R is the perfect gas constant, „ _ ^ 1 " with y the Euler 

constant (^=0.577216) and M the atomic weight (gmol 1 ) of the 

atomic species M. In table 2, we give the values of the parameters 
A and B that we used in our calculations. For the other collisions 
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between neutral and charged particles, we considered them as 
elastic and that the charge particle engenders a dipole in the 
neutral particles during the collision. From the paper of Kihara 
and al, we obtain the cross section by the formula [14,18, 23, 25]: 
4(/ + l) 
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Av are parameters given in [18]. The polarisabilities a of 

neutral particle that we use in our calculation are given in [11, 12, 
14, 19]. 

Table 2: Parameters A and B characterized resonant charge 
exchange used in the present work. 
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The Debye length taken in to account the presence of electrons 
and ions is given by [26]: 

(18) 
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Where £o is the vacuum permittivity and q; the electrical charge of 

the i th species. 

4) Results and discussions 

4.1) Test of calculation program 

In this section, we have tested our calculation program by 
comparing our results with a well-know values of dry air of the 
work of P. Fauchais and al [27]. On figures 1, 2 and 3, we 
respectively compare the total thermal conductivity, the electrical 
conductivity and the dynamical viscosity of our work with those 
of P. Fauchais. Figure 1 shows that our results of total thermal 
conductivity and those of P. Fauchais and al are in good 
agreement with a less than 20% maximum variation around 
4500K, 7500K and 9000K. The electrical conductivity in figure 2 
is in good agreement in all the range of temperature. For 
dynamical viscosity, the maximum discrepancy around 10700K is 
less than 10%. The different discrepancies in total thermal 
conductivity and dynamical viscosity can be due to the results of 
equilibrium composition and the integral collisions used in the 
transport coefficients calculation. 



Comparison of total thermal conductivity data of 
dry air plasma at atmospheric pressure 



t V 
/ \ 



v 



Species 


A 


B 


References 


H H + 


28.69 


1.3 


[12] [14] 


o o + 


19.5 


0.832 


[12] [14] 


N N + 


26.61 


1.27 


[12] [14] 


o 2 o 2 + 


24.05 


1.132 


[14] 


NO NO + 


23.61 


1.095 


[14] 


N 2 N 2 + 


24.5 


1.032 


[14] 
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Temperature (K) 

Figure 1: Comparison of thermal conductivity of our work with the one of P. 
Fauchais and al. 

Comparison of electrical conductivity of 
dry air plasma at atmospheric pressure 



3.4) Charged - charged interactions 

The average cross section for interactions between charged 
species has been calculated from coulomb potential [1, 5, 12, 14, 
25]. Using the expression from transport sections efficacies of 
Liboff, Devoto [1,5] established the following formulas: 



P. Fauchais 



/ 



/ 
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Figure 2: Comparison of electrical conductivity of our work with the one of P. 
Fauchais and al. 

Comparison of dynamical viscosity data of 
dry air plasma at atmospheric pressure 



P. Fauchais 
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Temperature (K) 

Figure 3: Comparison of dynamical viscosity of our work with the one of P. 
Fauchais and al. 

4.2) Influence of water vapor's rate on transport coefficients 

In figure 4, we have plotted the total thermal conductivity for 
different proportions of water vapor. Each curve exhibits two 
peaks: the first one around 3500 K is related to O2 and H 2 0 
dissociation, the second one around 7500 K to N 2 dissociation. 
The peaks related to ionization reactions occur at temperatures 
above 12000 K [1], which is out of our temperature range. The 
peak which appears around 3500 K, increases when the rate of 
water vapor increases in the mixture and the one around 7500 K 
decreases when the rate of water vapor increases in the mixture 
because the increase of water vapor's rate in the mixture 
contribute to the increase of the density of chemical species which 
dissociates at 3500 K (H 2 0) and decreases the chemical species 
which dissociates around 7500 K (N 2 ). The peak around 3500 K 
is very sensitive to the variation of the rate of water vapor. For 
temperatures below 6000 K and above 8000 K, the total thermal 
conductivity increases with the rate of water vapor in the plasma, 
and for temperatures between 6000 K and 8000 K, the total 
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thermal conductivity decreases when the rate of water vapor 
increases in the medium. 

The researches undertaken by various authors Hertz [1, 5], 
Vacquie [1,5] have showed that thermal conductivity has a very 
significant role on the time constant of arc (t). Consequently the 
shape of thermal conductivity curve gives important information 
on the performance of the gas used for the electrical arc extinction 
[1]. From the study of several gases, Hertz [1, 5] proposes an 
empirical relation linking the arc time-constant (x) to the 
maximum value of the reaction thermal conductivity component 
[1, 5, 28] of the gas: r.A™* = C' e . This relation shows the 

role of the reaction thermal conductivity component in the 
velocity of electrical arc extinction. The more the reaction thermal 
conductivity is high, the more the electrical arc extinct rapidly. 
So, the shape of the curve shows that the increase of water vapor 
in the mixture can be interesting for electrical arc extinction. 

In figure 5, we have plotted the electrical conductivity versus 
temperature for different proportions of water vapor. The 
electrical conductivity increases with temperature because the 
density of electrons increases when the temperature increases in 
the medium. The electrical conductivity decreases when the rate 
of water vapor increases in the mixture, because the increase of 
the rate of water in the mixture decreases the density of chemical 
species NO which are responsible by ionization of the production 
of electrons in low temperatures. That can be interesting for a 
good performance of the circuit breaker which uses this kind of 
gas. The shape of the curve shows that the electrical conductivity 
increases with the rate of water vapor in high temperature [1]. 

The dynamical viscosity of the studied plasmas is given versus 
temperature on figure 6. We obtain the shape of curves found in 
the literature. The results show an increase followed by a decrease 
versus temperature. The dynamical viscosity increases with 
temperature and reached its maximum value around 10000 K, 
which characterizes the transition between neutral or partially 
ionized plasma towards totally ionized plasma [13]. The decrease 
is due to the apparition of electrons and ionized chemical species 
[14] in the plasma. The dynamical viscosity decreases when the 
rate of water vapor increase in the mixture, because the molar 
weight of water molecule (H 2 0) is inferior to that of nitrogen 
molecule (N 2 ). We noticed a rapid increasing of the dynamical 
viscosity around 3500 K which correspond to the temperature of 
dissociation of water molecule (H 2 0). 

Thermal conductivity curves for P=1 atm 



20%Air-80%Water vapor 

— - 50%Air-50% Water vapor 

90%Air-1 0% Water vapor 



Electrical conductivity curves for P=1 atm 
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Figure 4: The total thermal conductivity versus temperature for the four mixture 
plasmas. 
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Figure 5: The electrical conductivity versus temperature for the four mixture 
plasmas. 

Dynamical viscosity curves for P=1 atm 
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Figure 6: The dynamical viscosity versus temperature for the four mixture 
plasmas. 

5) Conclusion 

We have presented the results of the transport coefficients of 
the plasma formed of air-water vapor mixtures. This study 
completes the one of Kagone and al [1] to have a database on the 
range of temperatures from 500 K to 30000 K for the modeling of 
plasma formed by air- water vapor mixtures. The results of 
thermal conductivity, electrical conductivity and dynamical 
viscosity versus temperature are presented and discuss. For 
temperatures less than 6000 K and above 8000 K, the thermal 
conductivity increases with the rate of water vapor. The peak 
around 3500 K is very sensitive to the variation of the rate of 
water vapor. The electrical conductivity decreases when the rate 
of water vapor increases at low temperatures. The dynamical 
conductivity decreases when the rate of water vapor increases in 
the medium. According to the empirical formula [1, 5, 28] 

T.A, e aL , = C' e and the decrease of electrical conductivity, the 

increase of the rate of water vapor can be interesting for the 
improvement of the performance of circuit breaker using this kind 
of gas. 
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